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Abstract 

We highlight the properties of a simple model (contained in our recent work) of 
the quantum dynamics of a single point vortex interacting with the nodal fermionic 
quasiparticles of a (f-wave superconductor. We describe the renormalization of the 
vortex motion by the quasiparticles: at T = 0, the quasiparticles renormalize the 
vortex mass and introduce only a weak sub-Ohmic damping. Ohmic (or 'Bardeen- 
Stephen' damping) appears at T > 0, with the damping co-efficient vanishing ~ T 2 
with a universal prefactor. Conversely, quantum fluctuations of the vortex renor- 
malize the quasiparticle spectrum. A point vortex oscillating in a harmonic pinning 
potential has no zero-bias peak in the electronic local density of states (LDOS), but 
has small satellite features at an energy determined by the pinning potential. These 
are proposed as the origin of sub-gap LDOS peaks observed in scanning tunneling 
microscopic studies of the LDOS near a vortex. 
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1 The Model 



A more complete discussion of the quantum dynamics of a vortex and nodal 
quasiparticles in a two-dimensional <i-wave superconductor appears in our re- 
cent papers [1,2,3] which also describe the connections to microscopic theory. 
Here, we focus on a simple low-energy model, which captures all the essential 
features. Our model has no zero-bias peak in the electronic LDOS at the vor- 
tex center (see also [4]), in contrast to computations in the traditional BCS 
framework [5,6,7]. The zero-point quantum motion of the vortex in the pinning 
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potential leads to sub-gap satellite peaks in the LDOS. Finally, as discussed 
in earlier work [8] (and not reviewed here) a proper accounting of Aharanov- 
Bohm-like phases in the vortex motion leads to periodic spatial modulations 
in the LDOS (see also [9]). These features have the potential to explain the 
key characteristics of scanning tunneling microscopic (STM) observations of 
vortices in the cuprate superconductors [10,11,12]. 

The degrees of freedom of our model are (i) the low energy S — 1/2 fermionic 
quasiparticles in the vicinity of the nodes of a rf-wave superconductor, de- 
scribed by the second-quantized Nambu spinor ^(r), and (ii) a (point) vortex, 
described by its first-quantized position r v and canonically conjugate momen- 
tum p v . In the absence of vortex pinning (considered later), the Hamiltonian 
has the very simple form (h — 1): 

H= ]T Jd 2 r^(r)H D (r,r v )V(r) , (1) 

nodes 

where H D is the Dirac Hamiltonian in the presence of a gauge field a [4]: 



H d = v 4^ + a *( r > r v) J ° z + v A (J^ + a y (r, r v ) j a x . (2) 



Here x denotes any of the four nodal directions in momentum space and y 
the direction perpendicular to it, a x,y,z are Pauli matrices, Vf and i>a are Fermi 
and gap velocities respectively, and the gauge field a is given by 



z x (r - r v ) 



so that it corresponds to a 7r-flux centered at the instantaneous vortex position 
r v ; in the Lagrangian formulation, Eq. (3) would appear as a Chern-Simons 
term. An assumption is that the vortex core size can be neglected in the 
cuprates, especially since there are no localized quasiparticle states in the d- 
wave vortex cores. The Hamiltonian (2) is derived from the full Bogoliubov- 
de Gennes Hamiltonian in the limit of an extreme type-II superconductor by 
application of the Franz- Tesanovic unitary transformation [4]. 

We have not displayed the 'Doppler shift' terms in Hp: their influence is con- 
sidered elsewhere [1], and they do not qualitatively modify the results below. 
Note also that a bare vortex Hamiltonian has not been explicitly included: the 
needed terms are generated by integrating out the quasiparticles. 
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2 Vortex dynamics in rf-wave superconductors 



First, we consider the influence of the nodal quasiparticles on the vortex dy- 
namics. This has been previously studied in a semiclassical theory [13,14]. Our 
results differ from this semiclassical theory which, we suspect, does not prop- 
erly capture the "quantum critical" aspects of the Dirac fermion dynamics. 

Instead of specifying a vortex Hamiltonian H v , we formulate an imaginary- 
time path-integral for quasiparticles coupled to a single vortex at the variable 
position t* v (t), and integrate out the fermionic quasiparticle fields. The result 
is an effective vortex action S v that contains only the quasiparticle contri- 
bution to vortex dynamics. We express S v in the frequency domain as an 
expansion in terms of the small vortex displacement from a fictitious origin, 
and calculate only the terms up to the quadratic order. Then, on general 
symmetry grounds we can write: 



FhHItvHI 2 + F ± {u)iz(r* v {u) x (r») 



(4) 



where F\\(u) captures the "longitudinal" dynamics, such as vortex friction 
and inertia, while F±(u) captures the "transversal" dynamics, such as the 
quasiparticle contributions to the Magnus force. 

The idea behind expanding the vortex action in terms of small vortex dis- 
placements is to make use of the solutions for the quasiparticle eigenstates ip n 
and eigenvalues e n in presence of a static vortex. Various virtual transitions 
between these extended quasiparticle eigenstates, caused by the vortex mo- 
tion, is what gives rise to renormalization of the parameters that characterize 
vortex dynamics. Thus, the contribution of quasiparticles to the vortex action 
at finite temperatures T = 1/(3 is found to have the following form: 



1 n,n> P « Cn-Cn'- ^ 

x \r v (u))U nyn ,\ 2 , (5) 

where /(e) is the Fermi-Dirac distribution function. The main problem is to 
calculate the transition matrix elements U n y = (ip n \ V|Vv)> and carry out 
the summations over the quantum numbers n and n'. Using the microscopic 
model given by (1) and (2), we find the quasiparticle eigenstates in presence 
of a static vortex, and substitute them into (5). In comparison with (4) we 
obtain at small frequencies: 

2 

F\\(u) = -ri\uj\+A 1 uj 2 \n(\uj\) + -^— + A 2 \uj\ 3 . (6) 
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At zero temperature the nodal quasiparticle contribution to the vortex mass 
is of the order of an electron mass: 



where A is the high energy cutoff of the Dirac Hamiltonian, and there is no 
vortex friction (77 = 0, A\ — 0) apart from a universal sub-Ohmic damping 
{A.2{v/\/v{) 7^ 0). These key results can be understood by a simple scaling 
argument [1] — non-analytic terms can only arise from infrared singular terms, 
and by power-counting, the first infrared singularity arises only at order a; 3 . 
Ohmic friction arises either at finite temperatures: 



or in presence of perturbations that create a finite density of states at zero 
energy p(0) 7^ (the scale T above is replaced by a scale oc p(0), which could 
be, for example, the strength of disorder, or Zeeman splitting). In general, 
quasiparticles do not contribute transversal vortex dynamics (F± = 0), unless 
p(0) 7^ 0. When the Doppler shift is included in calculations, these conclusions 
remain qualitatively the same in terms of their dependence on T, p(0) and 
va/v{, but the numerical co-efficients change significantly. 

Our findings are different from the semiclassical results in two important ways. 
First, we find a finite renormalization of the vortex mass due to quasiparti- 
cles, while the semiclassical approach predicts a mass that diverges in small 
magnetic fields as H^ 1 / 2 . Second, we find that quasiparticles cannot give rise 
to vortex friction despite their gapless spectrum, unless they are helped by 
a finite temperature or perturbations such as disorder. In contrast, the semi- 
classical approach predicts not only that friction is possible in the limit of 
infinite quasiparticle scattering time, but also that quasiparticles can produce 
transversal forces on the vortex. We hope that future experiments will resolve 
these discrepancies. 



3 Quasiparticle spectra in the vicinity of a vortex 

Now we turn to the converse problem of the influence of the vortex motion on 
the quasiparticle spectrum. From the results in Section 3, we know that the 
effective action for the vortex acquires a mass m v , and that effects of damping 
are negligible. We are interested in the electronic spectrum near a localized 
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vortex, induced either by the repulsive interaction with other vortices, or by 
a pinning potential associated with disorder. In both cases, for small vortex 
displacement, it is reasonable to make a harmonic approximation, and so we 
use the Hamiltonian H — > H + H v where 

The effective vortex mass m v , and the harmonic trap frequency uj v will be 
treated as given parameters, which can be determined microscopically [1,15]. 

We will apply several simplifications that amount to removing all sources of 
anisotropy in our model: we will set Vf — v&, use a single vortex harmonic fre- 
quency, and completely neglect the Doppler shift. Such simplifications are, of 
course, not quantitatively justified in realistic circumstances, but allow gain- 
ing a deeper physical insight and nevertheless produce quasiparticle spectra 
that are remarkably similar to the ones observed in experiments. 

3.1 Perturbation theory 

We define the operators 6^ and (/i G {x, y} that raise and lower respectively 
the quantum numbers of the vortex in the harmonic trap. By inserting the 
resolution of unity in terms of the trapped vortex eigenmodes \n x ,n y ) we 
can write the Hamiltonian perturbatively as H — H + H 1 + • • • , where the 
unperturbed Hamiltonian includes effects of the vortex zero-point quantum 
fluctuations: 

H = tu v b% + J dVttVo* , (10) 

while the perturbation describes resonant scattering of quasiparticles from the 
fluctuating vortex: 

H 1 = J d 2 r^ (y&t + h.c.) * . (11) 
Here we have introduced the following matrix elements: 

V^o(r) = (0,0|^(r)|0,0) 

V x (r) = (l,0\H D (r)\0,0) (12) 
V^(r) = (0,l|# D (r)|0,0> . 

The remaining terms in the full Hamiltonian correspond to scattering events in 
which the trapped vortex undergoes two or more virtual transitions between 
its discrete eigenmodes; such scattering processes are also generated in the 
perturbation theory, and their physical effects can be qualitatively obtained 
from (10) and (11) alone. 
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All calculations are performed numerically in the basis that diagonalizes (10). 
After all the simplifications, the quasiparticle eigenfunctions of H are char- 
acterized by the following quantum numbers: "charge" q — ±1, angular mo- 
mentum / G Z, and radial wavevector k > 0. Their spectrum e = qk is gapless 
at the gap nodes, and there are no localized states in the vortex cores of pure 
d-wave superconductors. The perturbation (11) describes coupling of Dirac 
fermions to a single bosonic two-dimensional oscillator. 

The quasiparticle LDOS is obtained from the quasiparticle Green's function 
Gu 

Muhtoi^) expressed in the spinor representation: 



with Ti tk (r) being the "Fourier" weight that translates between the position 
and momentum representations. Perturbative expansion of the Green's func- 
tion yields an expansion of the LDOS that satisfies the following scaling form: 



The small parameter a is the ratio of the perturbation energy scale and the 
vortex harmonic frequency: 



3.2 The quasiparticle LDOS 

The resonant scattering of quasiparticles from the fluctuating vortex leads 
to interesting effects already at the one-loop level. The quasiparticle LDOS is 
plotted in the Figure 1 as a function of energy at gradually increasing distances 
from the vortex trap center. The resonant peak height is proportional to a 2 at 
the trap center, and disappears over a length-scale comparable to the extent of 
vortex zero-point oscillations (m v c<j v ) -1 / 2 . This peak naturally lies at a sub-gap 
energy whenever the length-scale set by the superconducting gap (comparable 
to the vortex core size) is smaller than the amplitude of vortex quantum 
oscillations. Despite the simplicity of our model, the calculated LDOS has 
the same qualitative features as the LDOS measured in the STM experiments 
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Fig. 1. Energy scans of the full LDOS at gradually increasing distances r from the 
vortex core. The plots are offset vertically for clarity, starting from r = at the 
bottom, and moving up with increments Ar = §.2ui~ l for a 2 = 1 and Ar = O.SSlo^ 1 
for a 2 = 0.3 (Aruj v » (5a)' 1 ). 

An important feature of the LDOS is the absence of a zero-energy peak ap- 
pearing in earlier computations [5,6,7]. This appears to be a consequence of 
the smallness of the vortex core [16]. Even in our model, if the vortex core is 
larger than the spatial extent of vortex zero-point quantum fluctuations, the 
zero-energy peak appears in the quasiparticle LDOS [2]. Therefore, the STM 
experiments hint that it may not be the detailed structure of the vortex core 
that is crucial for shaping the quasiparticle spectra near the vortices in the 
cuprates; a viable alternative are quantum fluctuations of vortices. 

The one-loop correction p 1 (e, r) of the LDOS is calculated from the Fock self- 
energy: 



where the wiggly line represents a vortex excitation, and the straight line 
a quasiparticle. pi(e, r) is shown in the Figure 2 in order to emphasize all 
of its features. In general, there is a discontinuity of LDOS at the energy 
e = u y . A major peak appears at the energy e oc au v when it does not 
interfere with the discontinuity. The additional secondary features of pi(e, r) 
are too small to be seen in the full quasiparticle LDOS. The main peak moves 
toward lower energies when the amplitude of vortex quantum fluctuations oc a 
grows, indicating a short-lived resonant bound state between the vortex and 
a quasiparticle. 

Future even more sensitive experiments may provide more ways to detect con- 
sequences of vortex quantum motion in d-wave superconductors. The spatial 
extent of the LDOS modulations allows measuring the energy scale (hw v m v v 2 ) 1 / 2 , 
while resolving the thermally blurred discontinuities of the LDOS would allow 
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Pi/Wv a 2 =l Pi/Oh a 2 =0.2 Pi/Wv a 2 =0.05 




Fig. 2. The one-loop correction p\ to LDOS at the vortex center as a function of 
energy. These plots show evolution of pi(e) as the small parameter a 2 changes. The 
magnitude of p\ scales as a 2 . 

direct measurement of hw v . Combined knowledge of both energy scales reveals 
the vortex mass m v and the strength of the vortex trapping force (inter-vortex 
interactions), which in turn can be related to other parameters, such as the 
superconducting gap and magnetic field. Effects of the moderately strong Mag- 
nus force merely reduce to quantitative modifications of these energy scales. 

In closing, we reiterate that the above model can also explain [8] periodic 
LDOS modulations observed in recent STM experiments [11,12]. Such modu- 
lations appear in the LDOS over the region of vortex motion, and consequently 
the experiments lead to an estimate of m v [15]. 



4 Conclusions 

Contribution of the nodal quasiparticles to the vortex mass is of the order 
of an electron mass. Vortex friction arises only at finite temperatures, or in 
presence of perturbations that create a finite density of states at zero energy, 
such as disorder. Similarly, quasiparticles whose density of states vanishes at 
zero energy do not contribute any transversal forces on the vortex. Smallness 
of vortex friction and mass at low temperatures have important implications 
for the flux-flow in the "normal state" of d-w&ve superconductors. 

The quasiparticle LDOS in the vicinity of a localized quantum vortex shows 
remarkable similarity to the experimental STM observations. Absence of a 
zero-energy LDOS peak is an indication of the small vortex core. Small sub-gap 
peaks appear in the LDOS as a result of resonant scattering of quasiparticles 
from the quantum vortex. Some secondary weak features also appear in the 
LDOS, including discontinuities at energies that reflect the discrete spectrum 
of the localized vortex. The future STM experiments might be sensitive enough 
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to detect these features, and allow measurements of the vortex mass and vortex 
trapping forces. 
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